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ABSTRACT 

Spectral molecular line profile observations of star-forming molecular clouds sometimes 
show distinct red asymmetric double-peaked molecular line profiles with weaker blue 
peaks and stronger red peaks. For some star-forming molecular clouds, such molec- 
ular transitions with red asymmetric line profiles and blue asymmetric line profiles 
(i.e. blue asymmetric double-peaked molecular line profiles with weaker red peaks and 
stronger blue peaks) may coexist in spatially resolved spectral observations, while for 
others, such molecular transitions with red asymmetric line profiles may completely 
dominate in spatially resolved spectral observations. Blue asymmetric line profiles are 
usually interpreted as signals of central core collapses, while red asymmetric line pro- 
files remain unexplained. In this paper, we advance a spherically symmetric self-similar 
hydrodynamic model framework for envelope expansions with core collapses (EECC) 
of a general polytropic molecular gas cloud under self-gravity. Based on such EECC 
hydrodynamic cloud models, we perform tracer molecular line profile calculations us- 
ing the publicly available RATRAN code for star-forming clouds with spectroscopic 
signatures of red asymmetric line profiles. The presence of red asymmetric line profiles 
from molecular cloud cores indicates that EECC processes are most likely an essential 
hydrodynamic process of star formation. With spatial distributions, we explore various 
profiles of molecular lines for several tracer molecules in different settings of EECC 
dynamic models with and without shocks. 
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1 INTRODUCTION 

Collapses, expansions, shocks and turbulence are several im- 
portant dynamic features of star formation processes oc- 
curring inside molecular clouds. Such molecular cloud dy- 
namic characteristics may be revealed by rich diagnostics 
and comprehensive analysis of molecular spectral emission 
line profiles. An extensively discussed theoretical framework 
of for ming low-mass stars is the 'inside-out collapse' sce- 
nario (|Shulll977l : IShu. Adams fc Lizand [l987l ) . This model 
describes an isothermal self-similar dynamic solution that 
has a collapsing core surrounded by a static envelope with 
an expanding boundary engulfing more and more mass 
into the collapsed region. By adopting empirically inferred 
temperature variations to replace the constant tempera- 
ture, this dynamic collapse structure may lead to double- 
peak molecular line profiles with blue peaks stronger than 
red peaks (i.e. blue profiles), as revealed by spectral line 
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observations of some molecular globules in early stages 



of star formation (e.g. Zhou et alHl993l : ISaito et al. I Il999l : 



iHogerheiide fc Sandell|[2000l ) . Collapse solution based on dy- 
namics of general polytropic gas sphere can also lead to blue 
profiles, yet with a temperature variation in volved in the dy- 
namic model in a self-consistent manner (IGao, Lou fc 
I2OO9I ). 

Clearly, a collapse model is not a full story of star for- 
mation, as there exist observational signature s that cannot 
be accounted for in the collapse scenario (e.g. Wilner et al.l 



200c : lBelloche et al ]|2002l :l van der Tak. Caselli fc CeccarellH 



20051 ). One important signature is the detected molecular 
emission lines with red asymmetry, i.e. optically thick emis- 
sion lines are red shifted relative to optically thin lines from 
the same source. A statistical survey of the observed opti- 
cally thick molecular lines show that a quarter up to ^ 30% 
of al l sources show red asymmetry (e.g. .Mardones et al.l 
19971 : lEvansI l2003l : iFuUer. Wilhams fc SridharanI |2005| ). 
Among sources of red asymmetry, double-peak molecular 
line profiles with red peaks stronge r than blue peaks (i.e. red 
profiles) are further identified (e.g. IPark. Lee fc Mversll2004l : 
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iTafalla et allbOOel : IVeiusamv et aLlbOOSl ). More specifically, 
for a number of known cloud sources, both red and blue 
asymmetries are observed towards the same transitions of 
the same molec ules but at different beam offsets away from 
the centre (e.g. iTafalla et ai]|l998l : [Matthews eFaLl |2006| ) . 
While blue profiles are recognized as the signature of ra- 
dial infalls or core collapses, what then do these red profiles 
imply? 

Earlier radiative transfer calculations with param- 
eterized flow structures indicate that rotation and ra- 
dial outflo ws in molecular clouds may produce red pro- 
files (e.g. lAdelson fc Leun3 Il988l ). Dynamics of bipolar 
outflows has b e en extensively studied in the p ast (e.g. 
Shu et al.lll99ll. Il994l: iFiege fc HenriksenI |l996al : lOstrikeij 



19971 : iMatzner fc McKeell 19991 : IShang et ahlbood ). and their 



spectroscopic signatures (e.g. m olecular li ne profiles and ra- 
dio maps) are also explored (e.g. lFiege fc Hc nrikscn 1996bi). 
Being widely observed, these m olecular bipolar outflows 
have velocities > 10 km s"^ (e.g. IWu et al.|[2005l : ISu et all 
which are too large to account for the widely ob- 
served red asymmetry signatures indicating a typical gas 
flow velocity of ~ 1 km s~^ (e.g . [Mardones et al. 1 1 19971 : 



iFuUer. Williams fc SridharanI |2005| ). Of course, for special 
cases when bipolar outflows are oriented close to the 
plane of the sky, molecular line splittings can be smaller. 
Line profiles (as well as millimetre continuum maps) in a 
molecular cloud with bipolar outflow show an asymmet- 
ric spatial distri bution according to the axis direction of 
the outflow (e.g. Di Francesco et al] I2001I : [Matthews et al.l 
I2OO6I : |j0rgensen et al.ll2007D . Red proflles caused by rota- 
tion around the central core have a systematic displace- 
ment in the veloc i ty of local standard of rest (Vlsr) 
l|Di Francesco et all l200ll . Redman et al. 2004). The spa- 
tial distribution of red profiles is also as ymmetric accord- 
ing to the direction of the rotation axis l|Park et al.l 1 19921 : 
IZhou|[l995l ). Though generally considered to be part of con- 
tributions to the line broadening, turbulence in clouds is 
sometimes invoked to form line proflles with red or blue 
asymmetries, especially those more violent ones in the outer 
layers of clouds (e.g. Ossc nko pf 2002; Lee fc Kim 2009). 
Contraction and expansion motions caused by large-scale 
thermal pulsations (with a typical oscillation period of ~ 
10^ yr) in starless cores can also be origins of asymmet- 
ric molecular line profiles from molecular clou ds in their 
early s tages of star formation, i.e. starles s cores (iLada et al.l 
20031: iRedinan. Keto fc Rawhngj I2OO6I : iKeto et al.1 I2OO6I : 



Aguti et al.ll2007^ 



Theoretical models of envelope expansion with core 
collapse (EECC) for star formation represent a sign i flcant 
recent deve lopme i it (ILou fc ShenI l20o4 IShen fc Loul |2004 
iLou fc Gadbood : IWang fc L ou 2008, Yu fc Lou 2005; Yu 
et al. 2006; Hu fc Lou 2008). Such EECC dynamic process 
might affect the initial mass function of stars (e.g. Nakano 
et al. 1995) and the environment of star formation (e 



iMatzner fc McKee|[2000l : lMoraghan. Smith fc Rosedl2008 
The crucial question is whether such global dynamic struc- 
tures actually exist in star-forming clouds with sensible 
observational diagnostics such as characteristic features in 
molecular proflles? The main thrust of this paper is to show 
that theoretical molecular line profiles based on the general 
polytropic EECC solutions with collapses, expansions and 
shocks are able to explain observations of red profiles and 



can further provide plausible inferences of the star-forming 
region. 

This paper is structured as follows. We first present 
possible self-similar dynamic structures of collapse and ex- 
pansion solutions fo r general polytropic molecular clouds 
(|Wang fc Lod |2008| ) in Section 2. Molecular spectral line 
profiles are then obtained from radiative transfer calcula- 
tions in Section 3. Exploration of the EECC cloud conditions 
that generate red asymmetric line proflles and our perspec- 
tive of model applications are also presented therein. We 
summarize and conclude in Section 4. 



2 COLLAPSE AND EXPANSION DYNAMICS 

2.1 General Polytropic Hydrodynamic Models 

To describe star-forming molecular clouds, we adopt the gen- 
eral polytropic self-similar model framework of lWang fc Lod 
(|2008li ) but without the random magnetic fleld. In spherical 
polar coordinates (r, 9, nonlinear hydrodynamic par- 
tial differential equations (PDEs) for spherically symmetric 
molecular cloud dynamics are 
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where mass density p, radial flow velocity u, thermal gas 
pressure p, and enclosed mass M depend on radius r and 
time t; G — 6.67 x 10~^ dyne cm^ g~^ is the gravitational 
constant and 7 is the polytropic index. Equations ([1} and ((Sjl 
are mass and radial momentum conservations, respectively; 
equation (|3} is another form of mass conservation. Equation 
([4]) is the conservation of speciflc entropy along streamlines, 
which implies a general polytropic state equation (EoS) p — 
K{r, t)p'^ with K{r, t) being a coefficient that varies with 
both time t and radius r in general. 

These nonlinear PDEs allow self-similar solutions and 
the pertinent self-similar transformation is given below 
7 1/2. 
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(7) 



(3n-2)G ' 47rG 

where x is the self-similar independent variable combining r 
and t in a special manner, and a{x), m{x) and v{x) are di- 
mensionless reduced mass density, enclosed mass and radial 
flow velocity, respectively. According to mass conservation 
equation (|3]), the reduced enclosed mass can be expressed as 
m{x) = a{x)x'^[nx — ^'(a;)]. Being an important thermody- 
namic variable, the gas temperature T is given by the ideal 
gas law 
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where fee, A* ^-nd mn are Boltzmann's constant, mean molec- 
ular weight and hydrogen mass, respectively. We shall adopt 
^ = 1 for typical star-forming clouds. For a finite dm(x)/dx 
as X ^ 0, equation ((Tjl gives a central mass accretion rate 



Mo = k^'^t^^^-^^mo/G 



(9) 



with mo being the central reduced point mass enclosed. For 
n = 1, the central mass accretion rate Mo remains constant; 
for n > 1 and n < 1, this Mo increases and decreases with in- 
creasing time, respectively. Indices 7, n and q are related by 
general polytropic EoS (g)) with q = 2(n-|-7-2)/(3n-2). The 
case n -I- 7 = 2 features an EoS for a conventional polytropic 
gas and 7 = n = 1 describes an isothermal gas. Self-similar 
transformation equations ((Sjl — (O make it possible to cast 
nonlinear PDEs HJ — (0]) into nonlinear ordinary differential 
equations (ODEs) in terms of x which can be solved numer- 
ically with analytical asymptotic conditions and by taking 
care of the sonic critical curve (see Wang & Lou 2008 for 
details). 

Analytically, we have a static equilibrium solution for 
spherical cloud, namely, a singular polytropic sphere (SPS) 



i; = , a 



-1 -l/(n-3nij/2) 



2(2 - n)(3n - 2) 



-2/n 



2(2 - n)(3n - 2) 



-l/(n-3nq/2) 



{3n-2)/n 



(10) 



(e.g. Lou & Hu 2009). This static SPS solution may be help- 
ful for speculating the origin of those self-similar collapse 
and expansion solutions (see below). 

In the limit of a; — > +00, we have asymptotic similarity 
solution to the leading orders, viz. 

-2/n 



Ax 



nA 



(3n - 2) 



-f 2(2 - n)n''-'A'-"+^"'"'^ 



(n-2)/n 



+Bx 



(n-l)/n 



(11) 



where A and B axe two constants of integration, referred 
to as the mass and velocity parameters, respectively. In the 
other limit of a; ^ 0^, the asymptotic central free-fall solu- 
tion is 



2mo 



(371 - 2)x 



1/2 



(3n — 2)mo 
2^3 



1/2 



(12) 



where m = mo is the reduced enclosed mass m{x) = 



for a point mass at the very centre. 



2.2 General Polytropic EECC Solutions 

We mainly focus on self-similar dynamic solutions for en- 
velope expansions with core collapse in the form of free- 
fall towards the cloud centre (i.e. EECC solutions; Lou & 
Shen 2004 and Shen & Lou 2004). Relevant parameters of 
five such selected self-similar solutions are listed in Table [TJ 
where a;inf is the outgoing boundary separating the collapse 
and expansion regions; mo is the dimensionless central re- 
duced point mass; and mtot is the total reduced enclosed 
mass of a spherical cloud with an expanding outer edge at 
a; = 5 for the star-forming cloud. We also list in Table [T]two 
important mass and velocity parameters A and B, which 



characterize asymptotic dynamic behaviours for x +00 
and serve as asymptotic 'boundary' conditions in construct- 
ing general polytropic EECC solutions (see eqs (26) and 
(27) of Wang & Lou 2008 for details). Solutions IV and V 
are two EECC shock solutions with 7 = 1.2 and n = 0.8, 
the upstream point x\ and downstream point X2 correspond 
to the same shock radius Rsh', v\ and 112 are the upstream 
and downstream reduced radial flow velocities, with nega- 
tive value being infall and positive value being expansion, 
respectively. 

Each of these selected EECC solutions has a central 
collapsed core and an outer expansion envelope with an 
outgoing interface at a;inf separating the two zones. At the 
beginning of evolution {t 0^), this boundary radius 
is approximately zero according to self-similar transforma- 
tion ([5|, which indicates that collapse begins from the very 
centre of the cloud. As time goes on, this boundary ra- 
dius increases with a constant x^^i at a variable speed of 
Ub = Xinf according to equation ([5]), which means 

that more and more mass are enclosed into the collapse re- 
gion (see expression (O). At a given time t, mass densities 
and gas temperatures of all these EECC solutions increase 
towards the cloud centre. Shock solutions involve disconti- 
nuities of flow velocities, densities and temperatures across 
shock radius Rsh- Their dynamic structures are shown in 
Figs. [T] and [5] by adopting estimated physical scalings de- 
scribed below. 

In constructing self-similar hydrodynamic solutions, we 
have presumed a priori the possible existence of such form 
of similarity solutions under plausible asymptotic conditions 
and actually derive them analytically and/or numerically 
from nonlinear hydrodynamic equations satisfying relevant 
physical constraints. Meanwhile, it is of considerable interest 
to figure out even qualitatively how a molecular cloud sys- 
tem can sensibly evolve into such a self-similar phase given 
a certain class of initial and boundary conditions. This is a 
challenge especially in view of the existence of several pos- 
sible asymptotic self-similar solutions and requires a deeper 
theoretical understanding (e.g. Wang & Lou 2008). At this 
stage, we tentatively offer speculations on possible scenar- 
ios leading to self-similar EECC dynamic evolution invoked 
in this paper for modelling the dynamics of certain star- 
forming molecular clouds. 

Similar to stellar oscillations widely studied observa- 
tionally and theoretically, molecular clouds in SPS equi- 
librium pop when somehow perturbed may give rise to 
acoustic puls ations on much larger spat ial and tempo- 
ral scales (e.g. ILada et a l. 2003; Redman, Keto fc Rawlingj 
l2006l : lKeto et al.ll2006l : lAguti et al.112007 '). For simplicity, we 
may envision purely radial acoustic pulsations with possi- 
ble radial nodes in spherical molecular clouds; for exam- 
ple, such radial pulsations might be induced or excited by 
a sufficiently massive companion or transient object. With 
idealizations, such acoustic pulsations might persist period- 
ically for a long time in molecular clouds. Realistically, such 
acoustic pulsations might be 'damped' in one or two 'pe- 
riods' due to radiative losses as well as nonlinear effects. 
Among various pulsation phases, it would be possible to 
have a phase characterized by core contractions with enve- 
lope expansions. With such 'initial' conditions in molecu- 
lar clouds, the nonlinear evolution may eventually lead to 
core collapse under the self-gravity while the envelope ex- 
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No. 


7 


n 


1 


a;inf 


mo 


mtot 


A 


B 






X2 


V2 


I 


1.1 


0.8 


-1/2 


1.50 


1.49 


5.0 


3.5 


2.2 










II 


1.2 


0.8 





1.77 


2.32 


10.0 


5.0 


2.4 










III 


1.2 


0.9 


2/7 


1.80 


3.74 


20.8 


10.0 


2.4 










IV 


1.2 


0.8 





0.51 


0.207 


4.6 


2.68 


0.47 


2.77 


0.01 


2.57 


1.54 


V 


1.2 


0.8 





0.10 


0.0254 


9.2 


7.72 


3.70 


1.70 


-0.30 


1.68 


0.58 



Table 1. Parameters of self-similar EECC dynamic solutions without and with shocks for five Models I— V labelled by 
roman numerals on the left most column. Among the three scaling indices 7, n and q, only two are independent and all 
three are related by general polytropic EoS ^ with q = 2{n + '^ — 2)/(3n — 2). Three parameters x\^f, mo and mtot are the 
dimensionless infall radius separating the inner collapse and outer expansion regions, the reduced central point mass, and 
the reduced total enclosed mass for a model cloud, respectively. Two coefficients A and B are mass and velocity parameters 
in asymptotic 'boundary' conditions in constructing corresponding self-similar solutions as x — > -|-oo. Parameters xi, X2, vi 
and V2 are the upstream and down stream locations and velocities for solutions across an outgoing shock front (i.e. Models 
IV and V). 



pands into the surrounding interstellar medium (ISM). It 
is emphasized that no pulsations are necessarily persistent 
in this scenario. We speculate that this might evolve into 
self-similar EECC dynamic phase advanced in this paper. 
Following this scenario, pulsations of molecular clouds with 
different phases may evolve nonlinearly into a variety of dy- 
namic states. 

The EECC solution is necessarily consistent with the as- 
pect of energy conservation. For most low-mass star forma- 
tions, the Kelvin-Helmholz time scale tKH ~ GMtot/{RL) 
is longer than the dynamic times cale (pG)~^^^, because o f 
lower radiative efficiency (see e.g. iMcKee fc OstrikeillioOTh . 
Here Mtot, R and L are the total mass, the outer radius 
and the luminosity of the cloud, respectively. Because of ra- 
diative inefficiency, most low-mass star forming clouds need 
an extra means to carry out gravitational energy that are 
released during the central accretion. Global envelope ex- 
pansion, which is a sensible dynamic solution as already 
shown, serves as the extra energy release. From the consid- 
eration of energy conservation, the EECC shock solutions 
may potentially offer valuable clues to the problem that 
accretion rates derived from observed luminosity are much 
smaller than those expected according to their dynamic evo- 
lution (i.e . the luminosi ty problem) in low-mass star forma- 
tions (e.g. lKenvon et al., 1990: McKee fc Ostriker.i2007l ). For 
the physical scenario of variable central mass accretion rate 
and thus variable luminosity for the formation of low-mass 
stars in molecular clouds, we provide a general polytropic 
model explain the 'luminosity problem' (Lou & Dong 2009 
in preparation). 

2.3 Physical Properties of Molecular Clouds 

The reduced dynamic variables should be converted to phys- 
ical variables as applied to realistic cloud systems. 

The ty pical infall ra dius of a molecular cloud is ~ 0.01 — 
0.03 pc (e.g. lMversl[2005l ). or ~ 10^ -lO" AU. As the reduced 
infall radius is at a; ~ 1, we may choose the length scale as 

k^^^t" ~ 4 X 10^ AU (13) 

in self-similar transformation ([111 — ((7|. The outer cloud ra- 
dius is set at i? ~ 2 x 10^ AU. The number density 
at t he infall radius is estimated by ~ lO'^ — 10^ c m~'^ 
(e.g. iHarvev. Wilner fc Mverj l2003l : lEvans et aP l2009l l. A 
reduced number density of unity implies 



(47rG/.tmHt^)"^ 9 X 10" cm"^ , (14) 

leading to an estimated dynamic timescale. With parameter 
scalings (|13p and (|14l) for clouds, physical variables of a cloud 
can be expressed as follows according to equations ([5} — ([9]), 



namely 

r = 4 X 10^ X AU , (15) 

u = 0.213 v{x) km s"^ , (16) 

TV = 9 X 10* a{x) cm"^ , (17) 

M = 0.204 m(a;)/(3n - 2) Mq , (18) 

r = 5.33 Q(x)^"^m(x)'' K , (19) 

Mo = 2.26 X 10"'' mo Mq yr"^ , (20) 



where N — p/ (/imu) is the particle number density. We note 
that the dynamic timescale of a cloud is also automatically 
fixed, i.e. td ~ 2.2 x 10^ yr according to scaling estimate (|14p . 
From scaling estimates (|13p and (|14p , the sound parameter k 
is estimated by k^^^ = 3.76 km s""'^ for solution III and by 
fc^/^ = 65.7 km s"'' * for solutions I, II, IV and V; while the 
upstream k^^^ jumps up to k^^^ = 70.8 km s~" * and k^^^ = 
66.5 km s"''* on the downstream side of shock solutions IV 
and V, respectively. 

With scaling expressions (|13p and (|14p . we estimate 
physical properties of star-forming clouds as exemplified by 
Models I through V tabulated in Table [1 Figures [1] and |1] 
illustrate the radial profiles for physical variables of these 
model clouds in two sets. 

For those Models without shocks (i.e. I, II and III), the 
infall radius i?inf are similar to each other, and their cen- 
tral core mass Mo and total mass Mtot are also comparable. 
Here, the infall radius Rini may be fairly small at the onset 
of a cloud core collapse and expands to encompass more gas 
and dust particles into the collapsed region as time goes on. 
We note in Fig. [1] that the overall number density and tem- 
perature values increase from Model I to HI; and the temper- 
ature of Model HI increases slightly outwards at large radii 
(~ 10" AU). We note in Table [2] a gradual increase of the 
two masses Mo and Mtot from Model I to Model III, caused 
by the increasing mas s parameter A as shown in Table [1] 
[see equation (26) in l|Wang fc Loul lioosl ')]. Our numerical 
explorations reveal that as n -I- 7 increases, the presence of 
EECC solutions calls for increasing values of A. Another 
variation is that the ratio of central mass point to the total 
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No. 


7 


n 


Mo 


Mtot 


Mo (Mq yr-1) 


i?inf (AU) 


(AU) 


ni 


"2 


I 


1.1 


0.8 


0.75 Mq 


2.52 M0 


3.38 X 10~'^ 


6.0 X 10^ 








II 


1.2 


0.8 


1.18 Mq 


5.10 Mq 


5.30 X 10~6 


7.1 X 103 








III 


1.2 


0.9 


1.91 Mq 


11.00 M0 


8.41 X 10-*^ 


7.2 X 10^ 








IV 


1.2 


0.8 


0.106 Mq 


2.34 Mq 


4.73 X 10"^ 


2.04 X 10^ 


11.1 X 103 


0.002 km s-i 


0.34 km s-l 


V 


1.2 


0.8 


0.012 Mq 


4.68 Mq 


5.77 X 10-* 


0.40 X 10^ 


6.8 X 10^ 


-0.07 km 


0.13 km 



Table 2. Physical parameters of cloud Models I— V, with 7 and n being their polytropic indices and scaling indices, respectively. Two 
dimensional masses Mo and Mtot are the central point mass and the total cloud mass inside _R = 2 X 10'* AU, while Mq denotes the central 
mass accretion rate. Being the infall radius, Rinf is the boundary between the infall and outflow regions. Parameter R^^ is the outgoing 
shock radius for shock solutions in Models IV and V, and u\ and U2 are the upstream and downstream radial flow velocities, with negative 
values denoting inflows. 




Figure 1. Radial profiles for physical variables of three EECC 
dynamic models I, II and III, with different polytropic indices 7 
and n yet without shocks. From top to bottom are: radial flow 
velocity (positive values for infall), number density and tempera- 
ture profiles. The abscissa is radius R in 10^ AU in a logarithmic 
scale. Dashed, solid and dash-dotted curves are the EECC dy- 
namic solutions for Model I with 7 = 1.1 and n = 0.8, Model 
II with 7 = 1.2 and n = 0.8, and Model III with 7 = 1.2 and 
n = 0.9, respectively. The dotted horizontal line in the top panel 
is for the zero velocity line u = km s~^. The infall radii for 
Models I, II and III are _Ri„f = 6.0 X 10^ AU, Ri^i = 7.1 X 10^ 
AU and -Rinf = 7.2 x lO'^ AU, respectively. Other parameters for 
these general polytropic self-similar solutions arc summarized in 
Tables [Hand El 



Figure 2. Radial profiles for physical variables of three EECC 
dynamic Models II, IV and V, with and without shocks for poly- 
tropic index 7 = 1.2 and scaling index n = 0.8. From top to 
bottom are: radial fiow velocity (positive values for infall), num- 
ber density and temperature profiles. The abscissa is radius R 
in 10"^ AU in a logarithmic scale. Solid, dashed and dash-dotted 
curves are the solutions for Model II without a shock. Model IV 
with a shock at ij^h = H-l x 10'^ AU, and Model V with a shock 
at = 6.8 X 10"^ AU, respectively. The horizontal dotted line 
in the top panel is for the zero velocity line u = km s~^. The 
infall radii for Models II, IV and V are R^^i = 7.1 X 10^ AU, 
Rinf = 2.04 X 10^ AU and Ri^f = 0.40 x 10^ AU, respectively 
Other parameters for these self-similar EECC solutions are sum- 
marized in Tables [1] and [2] 



mass Mo/Mtot decreases from ~29.9% of Model I to ~23.2% 
of Model II to ~17.4% of Model III. This phenomena may 
be explained by accounting for general polytropic SPS so- 
lution (|10|) . of which smaller 7 and n values lead to mass 
distributions with more central concentration. 

By comparing Models II, IV and V, we can infer the 
effects of shocks in self-similar EECC dynamic solutions. 
From Fig. (2] we see that central infall velocities for EECC 
solutions with shocks (i.e. Models IV and V) are greatly sup- 
pressed, and their infall radii are smaller as compared 
to Model II without shock (Table [2]). As a consequence, the 
central mass accretion rate Mq and the central mass Mq for 
Models IV and V are smaller, taking only 10% and 1% the 
values of Model II, respectively. The low values of mass ac- 
cretion rates and central masses make these shock solutions 
potentially applicable for the formation of brown dwarfs. We 
should also note from Table [2] that the total mass Aftot for 



these three models are comparable, which suggests the possi- 
bility that under certain conditions protostars with different 
masses may form from molecular clouds with similar masses. 
In other words, molecular clouds of comparable masses have 
additional freedoms to give rise to central protostars of dif- 
ferent masses. The very low masses and mass accretion rates 
of EECC shock solutions are caused by efficient envelope ex- 
pansions with small inner boundary radius (Table Inside 
the shock radius i?sh, the particle number density A*' and gas 
temperature T of Models IV and V are higher than or close 
to those of Model II (see Figure [21). 

From Table [5] we get the information that the cen- 
tral mass of a cloud (i.e. protostellar mass) is not tightly 
or directly related to the total cloud mass for differ- 
ent model solutions, especially for EECC dynamic solu- 
tions with shocks. This conclusion agrees with the re- 
sult that the star forming efficiency Mo/Mtot is mainly 
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determined by the comparison between mass accretion 
and outflow rates (e.g. bipolar outflows therein), and is 
weakly dependent o n the core masfl Mtot and num- 
ber density A'' (e.g. [Nakano. Hasegawa fc NormarJ 1 19951 : 
iMatzner fc McKedbOOd '). However, a constant ratio between 
mass outflow and accretion rates is always presumed for a 
certain multi-star-forming cloud, and the conclusion that 
initial mass function (IMF) of protostars should be closely 
connected with the cloud mass function (CMF) is u sually 
drawn (e.g. iHennebelle fc Chabrieij|2008l : lMversll2008l ). But 
there is no obvious reason that protostars forming in the 
same cloud should have the same mass outflow ratio over 
accretion. Based on our analysis and results, we strongly 
suggest that solutions with different ratios between outflow 
and accretion exist in the same molecular cloud, and IMF 
of protostars can differ from CMF significantly, especially 
for those protostars with extremely low masses (i.e. brown 
dwarfs) . 

3 SPECTRAL LINE PROFILE SIGNATURES 

We perform radiative transfer calculations for molecular 
spectral line profiles using the publicly ava ilable numer- 
ical code RATRAN (|Hogerheiide fc Sandelll [2000) under 
spherical symmetry. This RATRAN code has been bench- 
marked among seven other radiative transfer codes by Van 
Zadelhoff et al. (2002) against examples of star forming 
clouds like B335. Based on the Monte Carlo method (e.g. 
Hammerslev fc HandscombI 1 19641 : IShreidej 1 19661 : iBernesI 
I979 I). this RATRAN code deals with both radiative trans- 
fer and non-local thermal equilibrium (non-LTE) excitations 
of atomic and molecular lines, making it readily adaptable 
to realistic astrophysical cloud systemslf] We adopt dynamic 
cloud models for line profile calculations. The dynamic and 
thermal parameters, namely radial fiow velocity u, number 
density N and gas kinetic temperature T are obtained con- 
sistently from the general polytropic EECC dynamic Models 
I-V. The temperature of dusts in a cloud is assumed to fol- 
low the gas kinetic temperature T. From the cloud centre to 
its outer radius R = 20000 AU, a molecular cloud is divided 
into 12 shells with enough accuracy for calculations (8 and 
10 shells are also tested separately but with little variance in 
line profiles) . However, physical properties of these shells are 
not uniform in order to show the transition between infall 
and outflow for all the dynamic models. 

In Figs. |3] |4] and O we show a sample of computed 
molecular line profiles, viz. HCO+ J= 1 - at 89.19 GHz 
and HCO+ J= 3 - 2 at 267.56 GHz, for the five dynamic 
EECC solution Models I— V. Molecular line profiles for tran- 
sitions CO J= 2 - 1 at 230.54 GHz, C^^O J= 1 - at 
109.78 GHz, CS J= 2 - 1 at 97.98 GHz and N2H+ J= 1 - 
at 93.13 GHz are shown in Figs. |6] and [T] All molecular 
data are obtained from the Leiden Atomic and Molecular 

^ We use the two terms 'core mass' and 'cloud mass' interchange- 
ably in the loose sense. 

^ For the purpose of testing and checking, our simple model ra- 
diative transfer calculations (Gao, Lou & Wu 2009) are also per- 
formed by adjusting optical depths in parallel to the RATRAN 
code calculations for samples of red profiles. The relevant results 
are comparable. 



Database (|Schoier et al.l I2OO5I ). Emission line profiles with 
different impact parameters b [i.e. distance of line of sight 
(LOS) from the cloud centre] and under different conditions 
of micro turbulence are shown in these figures. More com- 
pressed separations between impact parameters for point- 
ings are chosen for LOS passing through the most inner 
regions of clouds, in order to clearly show the transitions 
between blue profiles and red profiles for Models I— HI. 
The number density of HCO"^ molecules is assumed to be 
proportional to the overall molecular number density with 
a constant ratio, i.e. A'jjco+ = 2 x 10~^A'^. Molecules in 
Figs. [S] and [7] are also assumed to be of constant abun- 
dance ratios, viz. Nco = 5 x 10"^iV, A^ciso = 1 x 10"'' A^, 
iVcs = 3 X lO-^iV and iVN2H+ = 1-5 x 10-^°N, respec- 
ti vely. All these molec ular abundances are chosen according 
to lTafalla et al.l (|2006l ). and the central abundance hole de- 
rived from fitting procedure therein is not assumed in our 
model calculations. Strictly speaking, for all molecules, the 
variation of abundance ratio should be seriously take n into 
account (jRawlings fc Yatesll200ll : iTsamis etallbOOSl ). and 
the constant ratio abundance adopted here is just a first 
order approximation. 



3.1 Red Profiles from EECC Dynamic Models 

The most distinct spectroscopic signature of EECC dynamic 
models with core collapses and envelope expansions is the 
double-peak molecular line profile with the red peak being 
stronger than the blue peak, referred to as the red profile. 
We readily see such red profiles from cloud models described 
by each EECC solution in Figs. [HHSl A more generic type of 
red profiles has no obvious central dip but shows a stronger 
red shoulder, which exist over a larger range of impact pa- 
rameter b for the LOS. In Figs. |5] and [T] red profiles also 
widely exist except for CO and C^**0 transitions in Model 
HI (Fig.EJ. 

We demonstrate two types of self-similar EECC dy- 
namic solutions according to the appearance of red profiles 
for molecular spectral lines from star-forming clouds. For 
Models I, II and HI, red profiles only exist for a fairly large 
impact parameter b, i.e. in Fig. |31 fe > 1.08 x 10^ AU for 
Models I and HI, and b > 1.35 x 10^ AU for Model II; while 
inside these radii (i.e. impact parameter), blue profiles for 
molecular lines emerge because of the more dominant role of 
core collapse towards the centre. This type of spatial tran- 
sitions from blue profiles to red profiles can be seen in some 
molecular clouds with spatially resolved observations (e.g. 
iTafalla et all I2OO0I : IWard- Thompson fc BucklevI I2OOII ). al- 
though not sufficiently spherically symmetric. For Models IV 
and V, whose infall radii are very small at the chosen epochs, 
red profiles for molecular lines exist throughout pixels with a 
range of impact para meter b; this can also be se e n in some of 
recent observations l|Thompson fc Whitell2004l : lAguti et al.l 
l2007f ). These results of our numerical exploration clearly in- 
dicate that the broad existence of red profiles is a character- 
istic signature of global envelope expansion, in contrast to 
the blue p rofiles which are characteristic feature of core col- 
lapse (e.g. IZhou et"al]|l993l : [Cao. Lou fc Wull2009l ). For all 
these EECC dynamic models without or with shocks, molec- 
ular emission line profiles decrease in the overall magnitudes 
as the LOS departs away from the core centre and gradu- 
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Figure 3. Computed molecular line profiles of HCO+(1-0) at 89.19 GHz (solid curves) and HCO+{3-2) at 267.56 GHz (dashed curves) 
from individual pixels of spatially resolved emission lines for five EECC dynamic models, viz. Models I-V. The ordinate is the brightness 
temperature in Kelvin and the abscissa is the projected velocity component along the LOS in the local standard of rest (LSR) in unit of 
km s~^. Rows from top to bottom show molecular line profiles for different dynamic Models I— V, respectively. Panels from left to right 
present LOS with seven different impact parameter h (i.e. 0.27 X 10^ AU, 0.54 X 10^ AU, 0.81 X 10^ AU, 1.08 X 10^ AU, 1.35 X 10^ AU, 
3.24 X 10^ AU, 5.94 X 10^ AU). Impact parameters for pointings are not uniformly chosen on purpose such that the transitions between 
blue profiles and red profiles can be clearly identified for Models I, H and HI; meanwhile, the spatial coverage is sufficiently large to 
show variations of molecular line profiles. Containing both turbulence and thermal effects, the intrinsic line broadening is chosen as 
A« = 0.17 km s~^ 



ally disappear as cloud density and temperature decrease 
further. 

By comparing molecular line profiles among Models I, II 
and III, we find slight differences caused by different EECC 
dynamic profiles. The higher spectral amplitude in Model 
III appears to be caused by the higher number density and 
gas kinetic temperature involved (Fig. [T|. The central dips 
for line profiles of Model III are shallower than those from 
Models I and II; this is probably because of the gradual 
temperature rise at outer radii of Model III (Fig. [TJ . How- 
ever, we cannot immediately claim that these differences in 
molecular line profiles are caused by different values of poly- 
tropic index 7 and scaling index n in dynamic models, as 
for the same set of indices, different asymptotic boundary 
conditions will also le ad to considerab l e variations in dy- 
namic profiles (s ee e.g. lLou fc GaolbOOdlWang fc Lou"2008l: 
iHu fc Lou|[200a '). For Models IV and V, molecular emission 
lines with red profiles present for all values of impact param- 
eter 6. Then how could we understand the role of expanding 
shocks in these self-similar EECC dynamic solutions? Dy- 
namically, in gas clouds with relatively smooth or less dras- 



tic mass density profile, shocks are more likely to happen, 
which push more mass outwards and lead to more effective 
envelope expansions. The broad existence of red profiles just 
represents small infall radius at the epoch and highly effi- 
cient envelope expansion. Therefore by referring to underly- 
ing dynamic models, the presence of red profiles for nearly 
all impact parameter h indicates that the central protostel- 
lar mass is very small (i.e. ~ 0.106 Mq in Model IV and 
~ 0.012 Mq in Model V) and the molecular cloud may form 
a brown dwarf at the centre. 

3.2 Effects of Optical Depth and Turbulence 

Proper optically thick conditions (i.e. absorption and scat- 
tering) in star-forming molecular clouds are responsible for 
such asymmetric spectroscopic signatures as in Figs. [SHU 
Both HCO^ J= 1 — and J= 3 — 2 transitions are deeply 
self-absorbed (e.g. iTafalla et alll2006l '). while the J= 3 - 2 
line transition has a lower optical depth and source function 
because of lower level populations on J=3 and J=2 at such 
a cold environment with T '--^ 10 K. This contrast causes 
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Figure 4. Similar RATRAN computations for molecular spectra line profiles for five EECC dynamic Models I— V as in Fig. El but with 
intrinsic line broadening Am = 0.5 km s~^. This intrinsic line broadening contains both turbulence and thermal effects. 



two effects on molecular line profiles: first, the intensities 
of J= 3 — 2 lines are weaker; and secondly, asymmetries of 
molecular line profiles are less apparent for J= 3 — 2 transi- 
tions. This comes to a widely tested conclusion tfiat optically 
thick transition lines offer diagnosis for the large-scale ther- 
mal and dynamic structures of molecular clouds. CO and 
CS transitions in Figs. [6] and [7] are also examples of opti- 
cally thick transitions, which show distinct red profiles. In 
contrast, C^**0 and N2H'*' transitions are examples of op- 
tically thin transitions, which simply show single peak line 
profiles. 



Turbulence in molecular clouds is another important as- 
pect that affects molecular lin e profiles in a significant man- 
ner (e.g. ^Arons & Max 1975: Larson '1981'; 'Lou fc Rosiie^ 
1 19861 : [Zweibel fc M cKcc 1995; MacLow 1999). In this pa- 
per, the effect of turbulence is subsumed in the form of the 
intrinsic line broadening Au, where another contribution to 
the broadening is thermal broadening of Aiithc ~ 0.1 kms~^ 
for clouds with a temperature of ~ 10 K. We show radiative 
transfer results for a smaller or comparable turbulent broad- 
ening in Figs. |3] and |4l and molecular line profiles for clouds 
under stronger turbulence are shown in Fig. [S] Besides the 
increase of line width, line profile asymmetry decreases as a 
result of the enhanced turbulence. 



3.3 Star-Forming Molecular Clouds 

By fitting spatially resolved spectral line profiles to real 
molecular clouds, we may infer their underlying dy- 
namic structures, and further estimate physical parame- 
ters (e.g. protostellar mass, central mass accretion rate, 
dynamic age etc.) of star-forming molecular cloud cores. 
There are several aspects for the data fitting: (1) selec- 
tion of self-similar EECC dynamic solutions, (2) adop- 
tion of proper physical scalings derived from empirical in- 
formation and (3) choice of suitable turbulent broaden- 
ing. Additionally, sub-millimeter continuum observations 
can serve as a constrain t on the rad i al profiles of den- 
sity and temperature (e.g. lAdamslll99ll : IShirlev et al.ll2000l : 
iHarvev. Wilner fc Mver7l2003h before fitting to molecu- 
lar spectral emission line profiles. For describing more 
realistic cloud situations, variations of molecular abun- 
dance ratio should be taken into account. Star -forming 
mol ecular clouds L1 517B (e.g. Tafalla ct al.l l200^ . L1544 
(e.g. Ivan der Tak. CaseUi fc CeccarellL2005i'l L1551NE (e.g. 
Moriartv-Schieven. Butner fc Wannierl Il995h . L483 (e.j 



Park et al.l |200GI : iTafalla et al.l |200G| . Carolan et al. 2008) 



are among the good candidates of molecular clouds likely 
involving EECC dynamic motions. 

In addition to red profiles for molecular lines, there are 
several other clues when applying EECC dynamic models to 
molecular clouds. Global cloud systems have an estimated 
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Figure 5. Similar RATRAN computations for molecular spectral line profiles for five EECC dynamic Models I— V as in Fig.O but with 
intrinsic line broadening An = 1.5 km s~^. This intrinsic line broadening contains both turbulence and thermal effects. 



speed range of 0.1 — 1 km for flows, which is smaller 
than bipolar outflows of typical speed > 10 km s^^. Ob- 
served with sufficient spatial resolutions, molecular line pro- 
files from clouds with EECC dynamics will maifest a circular 
symmetry; while bipolar outflows have line emissions with a 
bipolar asymmetry in spatial distributions, as we mentioned 
in Section 1. Frequency resolution of about 0.1 km s~^ is 
needed to resolve spectroscopic signatures as shown in Figs. 
[3] to [5] High spatial resolution (e.g. 2" or smaller for molec- 
ular clouds at ~ 200 pc) is also very important for detecting 
variations of molecular line profiles as a function of radius 
in cloud core, which helps to distinguish different underlying 
dynamic structures. 



4 SUMMARY AND CONCLUSIONS 

We invoke self-similar general polytropic EECC dynamic so- 
lutions without or with shocks to model the global evolution 
of a certain class of molecular clouds. By specifying relevant 
parameters plausibly estimated for molecular clouds, we il- 
lustrate several examples of general polytropic EECC cloud 
solutions. On the basis of these cloud solutions, we perform 
radiative transfer calculations to produce molecular line pro- 
files to confirm the viability of EECC model framework. 

Through extensive numerical explorations, we demon- 
strate that the widely observed 'red profiles' in molecular 



emission spectral lines from star-forming clouds may well 
serve as important diagnostics for revealing the underly- 
ing EECC self-similar hydrodynamics in molecular clouds 
(Lou & Shen 2004; Shen & Lou 2004; Lou & Gao 2006; 
Wang & Lou 2008). From the point of view of general poly- 
tropic hydrodynamics and radiative transfer, our explana- 
tion for the mystery of 'red profiles' in emission spectral 
lines appears natural and physically sensible. In particular, 
a molecular cloud characterized by an envelope expansion 
with a simultaneous central core collapse represents novel 
scenario. Based on EECC solutions with or without shocks, 
optically thick molecular emission lines from gas clouds can 
show red profiles for all impact parameter h of LOS from 
the cloud centres, or just outside a certain b value (with 
the inner region showing blue profiles). The optical depth 
of transition lines and turbulent broadening caused by mi- 
cro gas motions will affect the appearance of red profiles. 
Different from those of bipolar outflows, emission lines from 
clouds under EECC dynamics will grossly show circular spa- 
tial symmetry, and the flow speed is typically smaller (e.g. 
~ 0.1 — 1 km s^^). By fitting spectral emission line pro- 
files of certain dynamic models with observed emission lines 
from star-forming clouds, it is possible to resolve the dy- 
namic structures of these molecular clouds. These processes 
will also give rise to more physical parameters of forming 
protostars and its cloud environment. 
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Figure 6. Model calculations for spectral line profiles of CO(2 - 1) at 230.54 GHz (solid curves) and Ci*0(l - 0) at 109.78 GHz (dashed 
curves) from individual pixels of spatially resolved emission lines for five EECC dynamic Models I— V. The ordinate is the brightness 
temperature in Kelvin and the abscissa is the projected velocity component in the local standard of rest (LSR) in unit of km s~^. Rows 
from top to bottom show molecular line profiles for different dynamic Models I— V, respectively. Panels from left to right present LOS 
with different impact parameter b, which arc not uniformly chosen such that spatial transitions between blue profiles and red profiles are 
clearly identified for Models I, II and III. The intrinsic line broadening is Au = 0.5 km s~^. 
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